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ABSTRACT 


Th«  growth  of  o vapor  bubhlo  In  a suparheatad  liquid  la  controlled  by 
three  factora:  tlw  Inertla.of  the  liquid,  the  aurfaee  tenaion,  and  the  vapor 
preaaure.  Aa  babble  growa,  evaporation  takea  place  at  the  bubble  boun- 
dary, and  the  tempexat  ire  and  vapor  pressure  in  the  bubble  are  thereby 
decreased,  Tha  heat  lallew  ;.‘equ.rement  of  evaporatien,  however,  depends 
on  the  rata  of  bubble  growth,  ao  ihat  the  dynamic  problem  ia  linked  with  a heat 
difluaion  problem.  Since  the  heat  dtffueion  problem  haa  been  solved,  a 
quantitative  formulation  of  the  dynamic  problem  can  be  given,  A aolution  for 
the  radius  of  the  vapor  bubble  aa  a function  of  time  ia  obtaiMd  which  ia  valid 
for  sufficiently  largo  radius.  This  asymptotie  solution  covers  the  range  of 
physical  interest  since  the  badiua  at  which  it  becomes  valid  ia  near  tha  lower 
limit  of  experimental  observation.  It  shows  the  strong  eUzct  of  heat  diffusion 
on  the  rate  of  bubhlo  growth.  Comparison  of  the  predicted  radius -time  be- 
havior is  made  with  axperimental  observations  in  superheated  water,  and  very 
good  agrsomsnt  is  found. 


I,  INTRODUCTION 


Whsn  the  vapor  pressure  ia  a liquid  eaceads  tha  ambisat  prasnurs,  it 
becomes  possrole  for  a vapor  buoMs  to  grow  from  a small  '’nuclaua"  ia  Vhe 
Uqiiid,  Thi  i anclous  is  a region  cl  noallquld  phase  and  pr  ssvunably  coaiista 
of  a 333  or  vhpor  phase  stahiUaed  on  a solid  particle.  Thi  rate  of  gro  ./th  of 
a vapor  b'lbhlj,  oacs  formed.  Is  delsririlaeQ  by  the  surlacs  tsasioa,  ihs 
llg’ud  iairlia,  aaJ  Vm  dVcJsTeacs  bilv/2;  'he  preasurr  ’.vithia  ths  bubbli  and 
t>,3  ambimt,  or  pTas3u*’e,  la  Ihs  iailial  stages  of  the  buVnla  sapaa- 

sion,  whjn  the  fore  33  are  marly  in  equilibrium,  the  growth  is  slo'.v  but  it  is 
accileratid  v/Uh  increase  ia  bubble  sine  because  of  tha  reducllon  ia  surface 
tension.  When  the  rate  of  bubble  growth  becomes  appreciable,  kov/ever,  the 
temperature  and  hence  the  praasura  within  the  bubble  drop  and  the  rate  of 
growth  Is  decreased.  One  might,  therefore,  expect  a maximum  la  the  velo- 
city of  the  bubble  wall.  The  reduction  of  the  temperature  within  tho  . - bble  Is 
a consequence  of  the  latent  h‘_*at  requirement  of  the  e vapor a^icr.  w'nic. 


plae«  at  the  vapar-liquid  Interfaca  as  the  bubble  grows. 

For  the  quantitative  solution  of  the  problem,  some  simplifying  physi- 
cal assumi^ions  may  be  made.  The  example  considered  in  detail  here  is  that 
of  a vapor  bubble  growing  in  moderately  superheated  water,  and  the  arguments 
for  tlw  validity  of  the  general  assumptions  are  Justified  lor  this  ease. 

It  will  be  assumed  that  the  bubble  is  spherical  throughout  its  growth. 

This  assumption  is  rsasonabla  providsd  the  radial  acceleration  and  vslocify  are 
small,  for  with  a spherically  symmetric  external  pressure  field  the  spherical 
shape  will  than  he.  stable  under  the  actfoncoC^urface  tension.  ^ 

Excluded  from  consideratioa  ts'the  asynntmetrie.  buoyant  force  of 
gravity  udiieb  becomes  importas^  if  the  bubble  growth  is  followed  for  so  long  a 
time  that  a eignlfieaat  translattonal  velocity  is  acquired.  A translational  velo- 
city of  the  bubble  as  a whole  not  only  causee  a deformation  in*  shape  but  also 
increases  the  rate  of  heat  inflow  to  the  bubble  over  that  used  in  the  analysis 
presented  below.  In  water  superheated  by  about  IQ^C,  no  great  error  is  intro- 
duced by  the  buoyant  force  provided  Uis  bubble  growth  is  3^ot  followed  beyond  a 
radius  of  appxoaiznately  1 mm. 

For  a superheat  of  lO^C  in  water,  a vapor  bubble  grows  from  its 

» 

initial  microscopic  siae  to  a radius  of  1 nun  in  a 'time  of  the  order  of  10  milU- 
scc.  The  corresponding  average  radial  velocity  of  10  cm/sec  is  very  sinall 
compared  wits  the  velocity  of  soxvad  la  thi  U-'juid,  so  that  compreasibiliiy  efiecls 
may  be  entirely  neglected  in  the  mottoa  of  the  liquid  producsd  by  the  bubble 
growth.  It  may  also  be  ohov/n  that  viscous  effects  are  uaimportant,  T-'.e 
hydrodynamic  equation  of  motion,  which  contains  thl  effect  of  tie  liquid  iasriia, 
is  thua  gr  latly  simpliiisd.  The  bubble  wall  veloc^  y is  also  very  small  oomeareu 
with  the  velocity  of  sound  in  the  vejor  so  that  compressibiUiy  effects  are  un- 
important here  as  well.  One  has  the  further  siinplUi cation  In  the  vapor  region 
that  the  ya.;or  density  is  so  small  that  its  inertia  offsets  may  be  neglected.  It 
then  follo'.vj  from  Bernoulli’ s equation  applied  to  a fluid  of  neglijihle  density 
that  the  pressure  a'/ithia  the  vapor  region  may  be  tahen  as  uniform.  In  addition, 
since  the  acoustic  velocity  in  the  vapor  ie  so  large,  the  pressure  within  the 
vapor  folio'./3  practically  inetantanaously  its  value  at  the  bubble  '.vail.  When  the 
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^lociiy  of  tho  bubblt  wall  ia  aulflctisntly  alow,  tha  proaauro  of  Ihe  vapor  is 
given  by  tho  equilibrium  vapor  pressure  of  the  liquid*  That  this  ia  the  casa 
for  the  present  problem  may  be  aoen  as  follows*  The  mean  velocity  appro- 
priate for  the  rate  of  evaporation  from  a liquid  surface  Is  (B  T/2«M)^' ^ 
triiere  B is  the  gas  constant*  T the  absolute  temperature*  and  M is  tho 
molecular  weight.  This  velocity  must  ba  reduced  by  the  coefficient  for  evapo- 
ration which  has  the  value  0i04  for  water.  ^ The  critical  vslocity  for  a water 
surface  at  approximately  lOd^C  is  therefore  about  3 m/aee*  which  ia  appre- 
ciably  greater  than  the  radial  yelociiiea  encountered  hare*  so  that  the  praeinure 
deBciancy  from  the  equiUbrium  vapor  pressure  may  be  neglected. 

I* 

The  tempefature  of  the  vapor  would*  in  general*  vary  with  position 
itt  the  bubble  as  well  as  Mth  time.  The  approximatiea  will  be  made*  however* 
that  the  thermal  diffusivity  D*  of  the  vapor  is  so  largo  that  temperature 
gradients  are  negligible  within  the  bubble.  In  water  impor  the  cl^ractariatie 
diffuaion  length  (20't)^^^  is  0.24  cm  for  tslo'^eec.  For  tha  graateat 
et^arhaat  conaidared  hsra  (iu  lO^C)  tha 'bubbla  radius  at  this  tima  ia  approx- 
imately 0.1  cm*  which  ia  aigsiflcantly  amallar  than  the  cliaractariaUc  diiluaion 
length.  The  approximation  of  uniform  trmperatura  within  tha  vapor  improvaa 
with  decreaae  ia  superheat. 

In  aummary*  the  physical  model  upon  v/hich  the  calculaSxoae  are  baaed 
conaiala  of  a apherlcal  vapor  bubble  *.vbich  haa  vuxiform  temperature 
praaaure;  the  temperatuve  of  the  vapor  ia  that  of  the  liquid  at  bubble  '.vail, 
and  the  pressure  ia  the  equilibrium  vapor  preaaure  for  that  tempi  .'alure.  In 
addition,  the  cTfecta  of  viacoaity  and  compreaaibiliiy  are  neglect?  1 both  in  the 
vapor  and  in  the  liquid,- 

The  equation  oC  motion  for  the  raJiui  R of  the  bubble  in  a naa/i  jcoua, 
incompreaaible  liquid  ae  a function  of  time  ie^ 

Rii  + ^ 

where  p ia  the  liquid  density,  la  the  e;ct?rnal  preaaure  in  the  liquid  or 
the  preasurs  at  Infinity,  and  p(R)  is  the  pressure  in  the  liquid  at  the  bubble 

^G.  Wyllie,  Pror,  Ro/i  Soc,  (London)  A 1*^7  ^ 333  (19-19). 

L.->i^b,  Hydrodyn  ii!'«io3  (Do\-i;r  I^uh.  Now  Yor!c,  19-13); 

M.  S,  Plo  Slot,  J,  Ar-ri.  M-.-ch.  16,  Z77(19V}), 
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bound»ryi  A dot  denotes  diflesentietion  virlth  respect  to  time.  The  pretiure 
it  sfiven  In  terms  of  the  vapor  pressure  p^  within  the  bubble  by 


P(R)  * Pv<T)-2e/R 


^2) 


where  v is  the  surlaee  tension  constant  and  p^(T)  is  the  equilibrium  vapor 
pressure  for  the  temperature  T at  the  bubble  boundary*  In  the  foUowins* 
the  small  variations  of  e and  p with  tompsratura  are  nsslected.  It  is  son- 
veniont  to  introduce  a radius  R.  defined  by  the  relation 


I 


2e/R, 


C3) 


Physically*  is  the  effective  Initial  radius  of  the  bubble;  it 


in  which  is  the  temperature  of  superheat  of  the  liquid  at  a distance  from 
the  bubble* 

represents  an  extrapolation  of  the  free  s^erical  bubble  down  to  the  equilibrium 

radius  for  the  siven  inittal  conditions.  It  should  lie  noted  that  a bubble  at  rest 

with  radius  R^  is  in  unstable  equilibrium.  The  actual  nucleus  from  uddch 

the  babble  {provrs  is  net  necessarily  spherical  and  its  surface  ensrijy  may  br: 

appreciably  Isss  than  4»eRj^  j hov/syar,  the  nucleus  from  which  the  bubble 

S70W3  and  the  free  spherical  bubble  of  radius  R^  sure  both  in  unstable  equili* 

brium  at  the  temp»jpaturc  T^  a».d  external  pressure  Table  1 g^ves  a set 

of  valuss  of  R^  at  various  superheat  temperatures  in  water  for  a pressure  oi 

P of  oTi  i atxn, 
o 

la  terras  of  the  parameter  R^,  Sq..  (1)  may  be  written 


RR  + R® 


P„(T)  - pjT  ) + (2t/R^){1  - R„/R) 


(^) 


or 


t » 


* “ . 3 • > 

— — — (R  R"")  s 

2R^R  dt 


Py(T)-p./T^)  + (27/R^)(1  -Ho/R) 


(3) 


If  the  coolinj  siiact  of  evaporaiion  is  dl  jrs^ardsd,  so  that  Py(T)  = p^(T^) , Sq. 
(3)  may  fcj  intsjrated  to  jive 


R*  » + (WSpRoJ  [*  * <2x/pR)  [l  -r//R^] 


(6) 


This  «otution  will  ba  rcforrad  to  as  Uia  Rayleigh  solutlon«  For  R>>  R^» 
£q*  (6)  becomee 


^ {Pv«V-**o] 


if) 


which  is  « constant* 


Ths  actual  moticm  daviatss  marhadly  from  that  pradictad  by  tha 

a 

Raylaigh  solution  becausa  of  the  cooling  effect*  Tha  heat  Q triiieh  must 
be  suppUed  to  the  bubble  per  onlt  time  is 


Q • (Wl)t.  . 


(a) 


adwre  L is  the  latent  heat  of  evaporatton  per  unit  mass  and  p*  is  tl»  rapor 
density.  This  heat  is  supplied  by  condnction  firom  the  liquid  into  the  bubble 
so  that 


Q * 4*a*k(dT/&r^  , 


(9) 


where  k is  the  thermal  condnetivity  of  the  liquid  and  Csr/dT}^^  is  the  temper- 
ature gradient  in  the  liquid  at  the  bubble  boundary.  Hence, 


(ST/ar).  . Js-  -4  — (31^P')  • 

* 3h  IT  dt 


(10) 


Tha  amall  variaUona  ol  L>  and  3t  with  teropaxaiure  ars  11  may  be 

shav/n  that  tha  contrihuttaa  cZ  tha  term  proporlioaal  to  — ”ich  lejs 

than  that  from  the  term  nrnpoytional  to  dll  /di,  Sqviation  (10)  tu/  he  approt- 
imated,  therefore,  by  O 


in) 


This  apprordmaUan  is  plausible  ph/jlcsUy  siacj  Ihj  vslvime  lncr3a'3T3  by 
orders  of  ma3nituJc,  v/hUe  ths  vapor  densii/  chanjea  very  liMVe  he^:-iuje  of 
the  small  temperature  decrease  during  the  bubble  jro'jrlh. 

Equation  (11)  tojether  with  the  specification  of  the  temperature  at 


^Thc  solution  for  the  motior4  of  n bubbl'i  u '1-r  auiii  constant  pr.:;:;'u; 
COr>'llhiQn3  wa3  wivon  by  Riyl;i;^K,  i'hil.  ’-I,  94  (1917)  ir-'i  applied  by 

to  th-‘  ca  ja  of  I.  roU  ;pji;i-j  b’;b’)l„-. 


a 1 ."o 
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iitftnlty  datormiaos  tho  tartperalura  Iteld  in  tha  liquid*  The  aoluiian  for 
this  toa^larat^ra  problam  with  tha  moving  boundary  R(t)  haa  baan  found^ 
under  the  aaauniptioa  that  tha  drop  in  tempo rature  from  to  the  value  T 

at  the  bubble  boundary  takea  place  in  a layer  of  liquid  adlacant  to  the  bubble 
which  haa  email  thiekneea  compared  with  R(t)«  Thia  approximation  of  tha 
'*thin  thermal  boundary  layer*‘  ia  juatifiad  phyaically  boeauae  the  thermal  diffu- 
aivity  of  tha  liquid  ia  amaUU  Tha  approximate  expreaaton  for  the  temperature 
at  the  bubble  wall  ia*^ 


a 

The  temperature  T ia  tfaua  given  in  terms  of  R and  R by  £q.  (12)  so  that 
Py(T)  ia  also  specified  in  termc  of  these  variables.  Equations  (S)  and  (12) 
therefore  determine  the  dynamic  problem  of  the  bubble  growth. 


H.  SOLUTION  OT  THE  EQUATION  07  MOTION 


Sqnatioaa  (3)  and  (12)  are  connected  whan  the  equilibrium  pres- 

sure ia  speciliad  aa  a functlo;.  of  temperabjre.  For  superbeate  not  too  far 
above  iha  hollinj  temperature  of  the  liquid  at  the  external  pressure  5^, 
the  vapor  praasura  may  be  approximatad  by  a linear  function  af  tha  ta  np-jra- 
turaj  • • •, 


Py<T)  - Pq 
P 

For  T = T^,  Sq,  (13)  glvsa 


r 1 \ 


2<r/j>a.  = A(T„-V 

because  of  Eq.  (3). 


(li) 


^M«  S.  Pie 3 set  and  S.  A.  Zwick^  J.  Appl.  Phya.  93  (1932). 

y 

Ref,  6,  Eq,  (20).  The  error  in  Eq.  (20)  may  be  estimated  from  Eq. 

(31)  of  thii  reforenca.  For  the  problom  d:3cuj;ijrl  here  it  waj  found  to  bo 

lodj  than  I0>3  of  tho  Hifforence  T - T ac  nnv  ti.no, 

* 


Equation  (S)  m&y  now  ba  Hrritton  aa 


I 


,3^2, 


h f.  • 


O »o 


2R‘*E  dt 

It  ia  convenient  to  use  In  place  of  R the  dimens lonlo as  volume  ratio 

p . , 


(15) 


<l6a) 


and  in  place  of  t the  dirvenatonleaa  variable 


U a (o/R 


*)  J R*(y)  <V  . 


(»6b) 


wheM  the  constant  a is  defined  aa 

r 1*/* 
• •[*] 


(Ifc) 


and  has  the  diznanaiona  of  reciprocal  time.  Equation  (15)  then  io  transformed 
to  ® 


J-.U 

o 


(u-y)^ 


(17) 


in  which  p*  =dp/du  and  the  dimensionleso  parameter  p is  siven  hy 
D 


A 

- n 


(li>).. 


Th-3  ph/iical  quantities  of  intsr  jst  are  now  expressed  as  follo*.vs; 


t a 

A " 
R ■ 


T-T_ 


if 


dv 

7^ 


V 


>/j 


(«R  /3) 


*%* 


• *-'  |> 


I 


* 


2 


9*M  dv 
(»^  • 


(Ha) 

(1W>) 

(He) 

(Hd) 


Tb«  valttaft  of  tha  vodiiu  R for  which  expcrimeatol  observations  nro 
rssdUy  ohtsiasbis  ors  moeh  grsotsr  than  R^  so  that  tbs  asymptotic  solution 
of  Sq.  (n)  (p»  1)  is  of  tbs  grsatsst  physical  interest.  It  Is  of  importance 
however,  to  saumias  tbs  initial  stassa  of  the  bubble  growth  since  a qualitative 
daacri^iioa  of  Ibis  portion  of  the  problem  is  naesaaary  for  an  understanding 
of  t!:s  aaymptotic  solution.  It  will  bs  shown  that  tbs  solutioa  in  tbs  asymptotic 
range  is  not  alfactsd  by  tbs  datails  of  tbs  miibsmaUeal  modal  ussd  to  ds  scribe 
tba  behavior  of  tbs  * .Nbls  for  R near  R^  so  that  tbe  uncertainty  in  tbs  sur* 
face  energy  for  ^jnall  vslues  of  R causes  no  dilJicultiss  in  tbs  physical 
apalieatioas  tr  bs  mads  bare. 


A.  laical  Gror/tb  ol 

it  is  ovl  Ijsl  i7ova  Sg*  (17)  that  a bubbl?  ol  radius  R^(y*  1)  ial "Ull/ 
at  reel  (p’(d)a  0)  will  remain  at  rest  in  ua  jVible  equilibrium*  Tbe  pree  net 
ma^bssae Ileal  mgdjl  reguirie  some  change  in  tbe  entemal  conJilioae  ^ » 
bv.Yele  growth  will  he  i:-dMaijrl, 

Tbe  peril sular  mechanism  v/hich  *>7111  be  considered  here  Is  a rise  in 
the  temperaturs  of  the  liquid.  Suppose,  for  saample,  that  there  ie  a constant 
heat  source  of  strengih  a per  unit  volume  so  that  the  amount  ol  heal  a Ided 
per  unit  volume  t)  from  t^O  to  t^t  is 


Tj  » at  , (20) 

Then,  if  the  temperature  of  the  bulb  liquid  is  at  t = 0,  at  time  t it  will  be 


(2’) 


T *•  ""  «“  *•  «*v«  by  «pl.clng 

in  Bq,  (12)  by  la  pUc«  ot  Bq,  (19d)  one  has 
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The  extension  of  Eq.  (17|  to  include  the  heat  source  term  is  thus 
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For  valoo.  of  p otar  H»  foJlUl  »alo.  of  oafty,  ooa  may  write 
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wbrro  lb.  .uballtulioa  tor  I faUo-.,a  from  E,.  (I9a).  Egoalioa 
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(23)  now 
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(24) 


Tbt  'onitaal  T f.  ,xl„a«iy  „„,3„  ^ 

liir.  rto.of  1 C/«c  <uri  B^-iO"  cm.  y 1.  approsimauiy  lo‘‘.  ih, 
•maiu...  of  lb.  conjiaat  y topli..  that  th.  torc.d  gro-.ylb  oway  from  Ih- 
oarlaM.  .qalUbrlum  polot  p , 1.  p.  . o i.  ,.ry  .low  „„m  bubbl , radio, 
h^a.  iocr.a..d  .oHloiootly  tor  tha  .urfaca  tan.ion  to  be  partially  ■■„Uxed... 
Thx3  imUal  slow  growth  is  a d-ilay  p.rio.l  in  Uto  bubbl  : .growth,  for  th  : 


bUbItVi  okftrtgdM  v«fy  tltilo  until  thn  t«maj|itthg  ietmi  bli  Iho  rlftht  alda  o£ 

D||I  (2:1^}  bouom«  abpfQciabioi  Tt\«.biibbU  liroWth  i«  ih«n  «o  rapid  (hat  the 
change  ib  the  bulk  temperature  o([  |he  UilUid  la  insignificant,  and  the  term  in  Y 
may  be  neglected.  The  delay  period  effectively  gives  an  initial  value  of  H 
elightly  greater  than  from  which  the  important  growth  begins*  An  approxi< 
mate  solution  may  be  found  for  this  initial  period  of  forced  growth  from  the 
equilibrium  point  by  linearizing  £q,  (23*)  ; tha't  is,  by  neglecting  terms  of  the 
second  order  in  p - 1 and  ita  derivatives.  The  details  of  this  calculation  will 
not  bo  presented  here  ; but,  as  would  be  expected,  the  forced  growth  away  from 
the  equilibrium  hae  an  exponential  behavior:  . 
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f I * —if 1 
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(25) 


The  temperature  at  the  bubble  wall  is  given  approximately  by 


T/5#  T. 


rR^Vj  zTO-t?)  ^ 
L A J P(3?^+l) 


CO 


(2S) 


The  constant  § is  that  root  of  the  equation 

- I s 0 

for  -ryhich  is  positive. 


3,  Gro-vth  Sehtivior  for  R ;?  '>Sq 

Aftor  the  bvaVol  * gro'vtii  has  bean  initiated,  there  ia  a rapid  xivi  in  the 
velocity  R until  the  cooling  jfect  becorors  irnporlant.  The  buhhie  v/all  velo- 
city djcreaaes  continuoujl/  thereafter.  No  o'o3ervatioaa  on  bubble  gro^vth  in 
this  region  are  aa  yet  availsble  and  the  details  of  this  aaelyeis  will  be  omitted. 

Of  pr.'ijent  intenei  Is  the  asymptotic  phase  of  the  bubble  gro  .vth  determined  by 
E^.  (17)  ’.vhicli  is  characterised  by  the  llmilinj  effects  on  R of  the  heat  diixu- 
oion  from  the  liquid  to  the  vapor.  As  R increases,  the  temperature  at  the 
bubble  wall  de  '>30333  steadily,  but  it  cannot  fall  below  ; for,  if  the  vapor 
temperature  f ^elow  T^,  the  pressure  difference  Py  “ would  become 
nejativi  and  t subble  gro',/th  would  be  arrested  and  eventually  reversed.  Such 


or 


behavior  i«  excluded  on  physical  grounds.  It  therefore  follows  that  the  in- 
tegral on  the  right  hand  side  of  £q,  (17)i  which  is  proportional  to  the  tempera- 
ture drop  at  the.  bubble  wall,  must  approach  a limit  as  t or  u oo  • A 
further  physical  argument  determines  more  precisely  the  asymptotic  behavior 
of  this  integral:  The  left  hand  side  of  Cq,  (17)  represents  essentially  the 
acceleration  effects  of  Uie  bubble  growth  in  the  tiqpiid.  As  the  bobble  grows, 
this  acceleration  tends  toward  zero  because  of  the  cooling  effect*  It  therefore 
follows  that 


1 (» 


as  a ^ 00 


(27) 


This  asymptotic  relation  may  be  inverted  to  yieldP 


P(»)  ^ 


as  tt 


»co 


(28) 


If  this  result  is  substituted  in  £q.  (17),  it  may  be  verified  that  Eqs*  (17), 

(27)  and  (23)  are,  in  fact,  consistent. 

Equation  (28)  is  not  yet  useful  since  it  provides  no  means  of  matching 
the  indicated  asymptotic  solution  of  £q.  (17)  with  a solution  valid  for  small 
values  of  p.  The  possibility  of  matching  solutions  depends  on  the  possibility 
of  shifting  the  asymptotic  solution  in  t (or  in  u)  so  as  to  account  for  the 
delay  period  in  bubble  growth.  It  is  necessary  that  one  be  free  to  shift  the 
asymptotic  solution  since  the  duration  of  the  deby  period  depends  completely 
on  the  choice  of  the  heat  source  term  '.vMle  the  subsequent  behavior  of  the 
bubble  is  independent  of  this  term,  A means  for  malcing  an  arbitrary  time 
shift  ia  furni  shed  by  noting  that,  in  addition  to  the  asymptotic  solution  of 
£q,  (2D),  ZSg,  (17)  also  possesaej  the  solution  p(u)sl.  It  v/ill  therefore  be 
assigned  that  the  asymptotic  solution  is  described  by 

p(u)  al,  0<U^u,  , 


Jp’(v)  dv  _ 

..  ■ 


I I 1 d 

1 - — 173  - - — — IP  . 
P ' 6p’ 


du 


u >u, 


(29) 

'30) 


O ,1  /» 

Equation  (27)  is  multiplied  by  (z-u)  ' and  integrated  from 
U = 0 to  u = z to  give  £q.  (28). 


-u- 


From  Eq.  (19a),  ona  hat  corresponding  to  Eqs,  (29)  and  (30) 


•o  that  Uj/n  represents  the  duration  o£  the  delay  period  in  the  growth.  The 
time  delay  may  be  introduced  esq^Ucitly  in  the  asymptotic  solution  for  u > u ^ 
by  use  of  the  fact  that,  if  p(u)  is  a solution  of  Eq.  (30).  p(u-f  u^)  is  also  a 
solution  udiere  u^  is  a constant.  A consistent  scheme  for  continuing  the 
asymptotic  solution  may  then  be  found  by  taking  the  solution  to  be  of  the  form 


P(tt) 


p(u)  « 1 , 
b. 


0 <u^tt. 


(32a) 


1 + 
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"s  .H 


2n(u-u^)"| 

C“-“o>  J ' 


u > u 
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(32b) 
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to  3JV30  When  the  coefficients  b|^  have  been  determined,  the 

diiijjreTice  (uj  -n^)  i*  fi^ad  by  the  requirement  that  p(U|)^  1,  The  delay 
period  is  then  determined  by  the  choice  of  • 

.When  (32b)  is  substituted  in  the  integral  on  the  left  side  of  Eq.  (30) 
the  result  is 


r slhL^  ^ 1 1 

J„  (a-v/^  3?^r(3/i 
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2P(-2/6)  •'5 
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(33) 


Q r 1**^  n/6 

. ^Higher  terms  are  of  the  form  I <n(u-u  )j  /(u-u  ) , where 

n arc  integers,  *- 


and 


m 


10, 


Sv^'^  ths  appundj.-c  for  th.;  evaluation  of  the  integral. 
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By  2q.  (30),  this  cxpreosion  i«  also  asymptotic  to 


It  (34)  is  axpanded  by  (32b),  th«  coefficients  of  corresponding  powers  of 
(u-u^)  may  bo  equated  to  give  a set  of  successive  equations  for  the  para- 
meters bj.  bj.  . , ,b^  , At  each  step,  one  has  a linear  equation  for  the 
unknown  parameter.  A tabulation  of  these  parameters  for  various  superheat 
conditions  in  watar  Ibr  an  external  pressure  of  1 atm.  i,  giv,n  in  Table  n. 

The  leading  terms  in  the  asymptotic  solution  arc  ^ 

P - u*/2  [l+0(u“*/^)]  . 

*<«|  [l^0(u‘/‘j]  , 


Thus, 


T-T„ ^ [l+0(a‘/‘)]  . 

[>  + 011-1/2,]  , 

a*R  * 

• -IT  [i  + o(i-'/*,]  . 


T-T  'w  . 
o 


•*»o* 


Th-  vapor  demity  pi  n»ay  be  w^pr^oximatad  by  a liasar  ItuicHon  of 
t>.5  tempsr-itu;:®  ov^r  ths  smiU  rari^e  of  pra.^ent  coaesrn.  Co«  V3.y,jnn/ 


P'-^C0.«1.  . ai  t->eo  , 

ky  2q.  (3i),  Th?  voluxn?  of  t>.a  b".Vol^  on  the  other  ho^id.  b/.iavji  lih> 

const.  . [l  + C{t’*/2)]  aj  t-->co  . 


aj  t — > CO  , 


It  l3  evident  that  the  vapor  ^nsity  varies  slov/ly  in  the  asymptotic  ran'r- 
compared  with  the  bubble  v ime  so  that  the  neglect  of  dpi/Jt  compared 

with  dn  /dt  in  the  formul  An  of  the  boundary  condition  of  E q.  (ll)  is 
juatifieds  / 


1 
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m,  COMPARISON  EXPERIMENT  AND  CONCLUSION 


With  th«  values  of  the  constants  given  in  Table  II,  raditts*time  curves 
have  been  computed  for  water  at  I atm  for  superheats  in  the  range  of  102”C 
to  106°C.  These  curves  are  presented  in  Fig*  1*  Inasmuch  as  the  delay  period 
may  be  chosen  arbitrarily  so  far  as  these  asymptotic  solutions  are  concerned, 
the  time  scale  is  determined  only  within  an  arbitrary  constant  which  varies  from 
one  curve  to  another.  The  actual  spacing  of  the  curves  as  presented  was  chosen 
so  that  the  time  intercepts  at  R a 0^  04  cm  were  equally  spaced. 

Observations  have  been  made  by  Dergarabedian^^  on  the  growth  of 
vapor  bubbles  in  superheated  waterA  The  comparison  between  the  theoretical 
curves  and  the  observed  values  is  shown  in  Figs.  2,  3,  and  4.  The  theoreti« 
cal  curves  were  obtained  by  interpolation  from  the  values  graphed  in  Fig*  1* 

The  time  origins  for  both  the  theoretical  curves  and  the  experimental  points 
are  both  arbitrary  so  that  a time  translation  of  the  theoretical  curve  has  been 
made  in  each  case  to  give  the  best  fit.  The  agreement  is,  however,  seen  to 
be  very  good.  The  importance  of  the  cooling  effect  is  evident  from  Fig.  2 
where  the  Rayleigh  solution  (Eq.  (7))  is  also  shown. 

Equation  (3S)  gives  for  the  leading  term  in  the  asymptotic  velocity 


jn 


(3») 


Whila  the  number  ol  terms  given  in  Eq.  (33b)  is  auiiicient  for  M;,h  accuracy  in 
the  range  of  values  of  R used  in  Figs.  1-4,  the  leading  term  in  Eq.  (33b)  is 
in  error  by  a^.prccilmaJ'.ely  10  to  30'/i  for  the  smallest  values  of  R,  This  lead- 

i t 

ii'.j  ierm  is  omJicijnlly  simple,  fcowevsT,  that  a rough  physical  argument 
should  give  the  essential  variation  of  3q,  (37),  The  argument  is  as  foUov/s; 
At  a time  t at  which  the  bubble  radius  R is  much  greater  than  R^,  the 
dilierjnci  behveen  the  temperature  Ir.  the  liquid  at  the  bubble  wall  and  that  in 
the  liquid  at  a distance  is  only  slightly  less  than  T^  temperature 

drop  ta!;es  place  principally  in  a liquid  layer  around  the  bubble  of  approximate 
thiclcness  given  by  the  diXfueion  length  The  heat  flow  into  the  bubble 


11 


P.  DcrgaraoccUT'.i,  J,  AppU  Mcch. , in  proas. 
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par  unit  tlm«  l«  thorefoY«  given  roughly  by 

4lfR*  . (38) 

(Dir* 

The  heat  requirement  per  unit  time  of  evaporation,  on  the  other  hand  is 

Q « L ^ (^»R^p*}^4vI^RLp>  . (39) 

When  (39)  ie  equated  to  (33),  there  results 

. (40) 

Lp»(Dir- 

ediich  agrees  in  order  of  magnitode  with  the  leading  term  as  given  by  Eq.  (37)« 

Some  experiments  have  recently  been  performed  by  Oergarabediaa  on 
vapor  bubble  growth  in  pure  CCI^  « At  moderate  superheats,  ita  vapor 
pressure  carve  Is  approsimataly  parallel  to  that  of  water,  displaced  to  lower 
temperatures.  If  the  rates  of  bubble  growth  in  ^se  two  liquids  are  com* 
pared  at  tha  same  value  of  the  temperature  difTerence  (T^  - Tj^),  they  should 
thus  be  rotijhly  in  the  same  ratio  as  lc{L>  p*  for  the  two  liquids.  This 

constant  is  about  3, 3 greater  in  vnter  than  in  carbon  tetrachloride, 
DergarabetUan'a  ohservatioas  on  bubble  gro'vths  are  in  good  agrseirent  with 
this  valus. 


-li- 
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APPENDIX 

Evaluation  of  tha  Tomperaturo  Integral  (Eg, 

By  difCorentiatittg  p(u)  (Eq.  (32b)  ) and  •ubstitutias  into  the  temperature 
integral  there  ro suits,  after  a change  of  variable, 
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Conaidar  a typical  integral  appeariag  in 
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For  0<a  <1,  3l'j(i)>9,  tiiii  ia  liyi^erjaom'jlric  f’jacUon 


^ riiOlCiaL  - 1 

P(3  + 1/2)  3 


m X — > 0 


(c) 


(a) 


By  an  argumant  bajsd  on  ths  theory  of  analytic  continuation,  however,  it  mny 
readily  be  ahown  that  (c)  is  valid  provided  only  that  < l»  « l^0»  *lf  -2, 

and  hence  that  (d)  holds  for  0<x<l,  Re(s)>-1,  (r  s =0,  (d)  i3 


E.  T,  Cop  i 'n,  Intro  iuctian  tv  TK:o.  y cl  Fuccti  ): 
Viriablj,  (Oxfo:-  1 Tr  - n,  Lo.v’  i, 
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monnlnglustt,  but  gives  a eorroct  result  under  a limiting  procedure.  By 
dUlerontlating  (d)  with  respect  to  s at  s s - l/z,  one  readily  finds 
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Arf  • 2ir  + ( 2a  Jt  + 2)  as  x —♦  0 


With  the  aid  of  (d),  (c),  £q,  (a)  therefore  becomes 
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afior  a rzarransamnt,  several  terma  cc^nesUing*  The  P -functions  appear> 

k > 

inj  in  (d)  have  bsjn  evaluated  numerically  in  (f).  By  comparison  \7ith  the 
asymptotic  form  for  p(u)  of  Sq.(30)  the  last  term  in  (f)  may  be  written 


(u-u 


; •>3-^p(U|) 


Inasmuch  as  p(uj)-  1 by  assumption,  Eq.  (f)  reduces  to  Eq.  (33)  of  the 
text. 
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Comparison  of  thaoretical  bubble  radiua-tims  values  with 
e:tperimental  values  hi  water  at  1 atm  with  a superheat 
temperature  of  103, 1 C.  The  solid  curve  is  the  Rayleioh 
growth  curve  neglecting  the  cooling  effect;  the  dashed 
curve  Includes  this  effect  by  the  method  indicated  in  the  text. 


I 
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Fig*  3 - Comparison  of  theoretical  radius-time  values  with  three 
sets  of  experimental  values  obtained  in  water  superheated 
to  104.5  C at  an  external  pressure  of  1 atm. 


